Aharonov-Bohm effect in undoped graphene: Magnetotrans- 
port via evanescent waves 
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Abstract. - Using conformal mapping technique, compact and general analytic expressions for 
the effects of magnetic fluxes on conductance and Fano factor of undoped graphene nanoflakes in 
pseudodiffusive regime are derived. 



Massless Dirac fermions in two dimensions have amaz- 
ing property, the finite conductivity (of order of conduc- 
tance quantum e 2 /h) in the limit of zero charge carrier 
concentration and no disorder [1]. It can be understood 
in terms of transport via evanescent waves (zero modes of 
the Dirac operator) [2]. Despite we deal with an ideal situ- 
ation (no scattering by defects) some transport properties, 
such as statistics of shot noise, are reminiscent of those for 
disordered metals [3] which justifies a frequently used term 
"pseudodiffusive regime". There are many theoretical [2-7] 
and experimental [8,9] works studying this regime in a con- 
text of graphene. This situation is very special. For con- 
ventional electron gas in semiconductors, in the absence of 
disorder, the states with definite energy (eigenstates of the 
Hamiltonian) can be simultaneously the states with def- 
inite current (eigenstates of the current operator) and it 
is the disorder that results in non-conservation of the cur- 
rent and finite conductivity. Contrary, for Dirac fermions 
the current operator does not commute with the Hamilto- 
nian ("Zitterbewegung") which can be considered as a kind 
of intrinsic disorder [2,10]. More detailed understanding 
of the pseudodiffusive transport in graphene is therefore 
of a general interest for quantum statistical physics and 
physical kinetics. 

Here we discuss peculiarities of the Aharonov-Bohm 
(AB) effect [11, 12] in the pseudodiffusive regime for 
graphene. The AB effect in graphene has been studied 
already theoretically [13-15] and experimentally [16,17]. 
Combining a general consideration of zero-energy states 
for massless Dirac fermions [18] with recent elegant treat- 
ment of the transport via evanescent waves by conformal 



transformation technique [6] we derive simple and gen- 
eral rigorous formulas for any graphene flake topologi- 
cally equivalent to the ring, avoiding both numerical sim- 
ulations and explicit solutions of Shrodinger equation for 
some particular cases. 

For the case of zero energy E = (undoped graphene) 
the Dirac equation 



a (-zV-A) i/j = Eip 



(1) 



for the two-component spinor ip, A is the vector potential 
and we use the units h = e = 1, splits into two independent 
equations for the components ip a (a = ±): 
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We will use the gauge VA = and introduce scalar mag- 
netic potential tp by 



thus, 



A -*tt A --^ 
Ax ~ dy' v ~ dx' 



V> = -B 



(3) 



(4) 



where B is the magnetic induction. The vector potential 
can be eliminated by a substitution [18] 
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Eq.Q being satisfied for an arbitrary analytic function 
X+ (z = x + iy) and complex conjugated analytic function 
X (z) ■ The latter can be found from boundary conditions. 
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To consider the pseudodiffusive transport one can assume 
that the leads 1 and 2 attached to the undoped graphene 
flake are made from heavily doped graphene, such that 
the Fermi wavelength of electrons within the leads is much 
smaller than any geometric size of the system under con- 
sideration. Thus, the boundary conditions read [2,3,5] 



1 



l-r = xp^ 

(2) 
+ 

(2) 



t = 1p 

t = v 



(6) 



where r and t are reflection and transmission coefficients, 
respectively, and superscripts 1 and 2 label the bound- 
aries attached to the corresponding leads. If the boundary 
of the sample is simply connected one can always choose 
tp = at the boundary and, thus, the magnetic field dis- 
appears from Eq.©; this fact was used in Ref. [7] as a 
very elegant way to prove that a random vector poten- 
tial does not effect on the value of minimal conductivity. 
Further we will consider a sample with a topology of the 
ring where the scalar potential tp is still constant at each 
boundary but these constants, tp\ and tp?. are different. 
Also, by symmetry, 

' W /x? ) =X<- 1) /x? ) 



(7) 



Further derivation follows Ref. [6] and the answer for the 
transmission coefficient T — |i| 2 has the form: 

Tj = 5 (8) 

3 cosh [2 (j + a) In A] V ; 

where j = ±1/2, ±3/2, ... labels zero modes of the Dirac 
equation, A is determined by a conformal transformation 
of our flake to the rectangle and 



P2 ~ Pi 

21nA 



(9) 



For simplicity, we will consider further the case of the 
Corbino disc with inner radius R\ and outer radius R 2 , 
when [6] 

A = y/R 2 /Ri (10) 

The conductance G (per spin per valley) and Fano factor 
of the shot noise F are expressed via the transmission 
coefficients (TSl) as 
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To calculate the sums in Eq. lfTT]) one can use the Poisson 
summation formula 



J2 exp[-2iirkx]f k (12) 



where n, k are all integers and fk is the Fourier trans- 
form of the function / (x) . Substituting Eq.® into (fill) 
one finds a compact and general answer for the effect of 
magnetic field on the transport characteristics: 



G 



F = 1— 



where 
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l + 2^(-l) fe cos(27rfca) 



fe=i 



(13) 
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(15) 



sinh (7r 2 fc/21nA) 

Eq.|4]) can be solved explicitly for radially symmetric 
distributions of the magnetic field B (r): 

R 2 r 

tp 2 - tp x = ± In (^) + J ±- J dr'r'B (/) (16) 
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where <I> is the magnetic flux through the inner ring. In the 
case of AB effect where all magnetic flux is concentrated 
within the inner ring one has 
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(17) 
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where we have restored world constants. 

Due to the large factor ir 2 in the argument of sinh in 
Eq. lfT5|) only the terms with k = 1 should be kept in 
Eqs. (fl3| and (fl4|) for all realistic shapes, thus, 
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where G = Kh^JW) - 

Oscillating contributions to G and F are exponentially 
small for very thin rings but, for sure, measurable if the 
ring is thick enough. For R 2 jR\ = 5 their amplitudes are 
5.3% and 40%, respectively. 

Pseudomagnetic fields describing by the vector potential 
can be created by deformations of the graphene flake [19, 
20]. Expressions lfl3"ll . lfl4|) can be applied also in this 
situation. 

Consider now a generic case with the magnetic (or pseu- 
domagnetic [19]) field B = within the flake. Then, the 
solution of Eq.Q can be obtained from the solution for 
the disc by the same conformal transformation which de- 
termines the function A. One can see immediately that 
Eq.lfl7|) remains the same. As for the expressions fl8|) . 
(191) they can be rewritten in terms of experimentally mea- 
surable quantity Go, 
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r-i + W-W*"® (21) 

where (3 — 2e 2 /hGo and we assume (5 -C 7r 2 . 

To conclude, conformal transformations [5,6] is a power- 
ful tool to consider pseudodiffusive transport in undoped 
graphene flakes of arbitrary shape, not only without mag- 
netic field but also in the presence of magnetic fluxes in the 
system. Experimental study of the Aharonov-Bohm oscil- 
lations and comparison with simple expressions ([20)1 , (|2Tj) 
derived here would be a suitable way to check whether the 
ballistic regime is reached on not in a given experimental 
situation. 
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